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We consider an extension of the holographic Ricci dark energy model by intro- 
ducing an interaction between dark energy and matter. In this model, the dark 
energy density is given by p\ = — ^aMpTZ, where 7Z is the Ricci scalar curvature, 
M p is the reduced Planck mass, and a is a dimensionless parameter. The interac- 
tion rate is given by Q = ^Hp^, where H is the Hubble expansion rate, and 7 is a 
dimensionless parameter. We investigate current observational constraints on this 
model by applying the type la supernovae, the baryon acoustic oscillation and the 
CMB anisotropy data. It is shown that a nonvanishing interaction rate is favored 
by the observations. The best fit values are a = 0.45 ± 0.03 and 7 = 0.15 ± 0.03 for 
the present dark energy density parameter Oao = 0.73 ± 0.03. 
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o\ '. 1 Introduction 
O 

Recent observations of type la supernovae (SNIa) have revealed that the present Universe 
is undergoing accelerated expansion pQ. This indicates that the energy density of the 
Universe at present epoch is dominated by dark energy with equation of state w\ = ^ < 
— |, where pa is the energy density of dark energy, and p\ is its pressure. Combining 
with observations of CMB anisotropy [2] and the baryon acoustic oscillation (BAO) [3J, 
the density parameter is determined as ^a = Pao/Pco ~ 0.7, where p cQ = 8.0992/i 2 x 
10 _47 GeV 4 is the critical density today and H = 100/i km/s/Mpc is the present Hubble 
parameter. 

The simplest way to explain dark energy is to introduce a cosmological constant jl]. 
However, this model suffers from two problems [5]. The first one is the fine-tuning prob- 
lem: the observed cosmological constant is extremely small compared to the fundamental 
Planck scale pa ~ 10 _120 Mp, requiring an incredible fine-tuning. The second one is the 
cosmic coincidence problem: why the cosmological constant and matter have comparable 
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energy density today even though their time evolution is so different. A variety of mod- 
els have been proposed to solve these problems, such as quintessence j6], phantom [TJ, 
quintom [S] and holographic dark energy [HI QUI Ell E2] ■ 

In particular, the holographic dark energy (HDE) models have been discussed exten- 
sively in recent years. These models are motivated by the holographic principle of quantum 
gravity [13]. From the condition that a system with size L would not form a black hole, 
it is required that the total vacuum energy should not exceed the mass of the black hole 
of the same size. Therefore, the dark energy density pa must satisfy L 3 p\ < LM 2 [T4"] . 
where M p = 1/ V^ttG is the reduced Planck mass. Saturating this inequality, pa is defined 

i'.v M 

p A = 3c 2 M p 2 L" 2 , (1) 

where c is a dimensionless parameter. In Ref. [13], the size L, which is regarded as an 
IR cut-off, was chosen to be the inverse Hubble expansion rate L ~ H^ 1 to naturally 
explain the observed vacuum energy density pao ~ "iM^H^. However, it was shown that 
this choice can not explain the accelerated expansion of the Universe at present [15] . This 
problem was solved in Ref. [9] by choosing L to be the future event horizon Rh . Extensive 
studies on this model have been done, and the HDE model with the future event horizon 
as the IR cut-off is found to be consistent with current observational data |16j . Various 
aspects on the HDE model have been discussed in Ref. [17] . There have also been further 
developments on the HDE model by introducing an interaction between dark energy 
and matter [JSJ [19]. However, it was pointed out that this model with L = Rh has a 
conceptual problem that the future event horizon, which determines pa, depends on the 
future evolution of the Universe, hence violates causality [20J. 

Subsequently, inspired by the HDE model, the holographic Ricci dark energy (RDE) 
model [10] was proposed in which pa is proportional to the Ricci scalar curvature 71. It 
was shown that this model does not only avoid the causality problem and is phenomeno- 
logically viable, but also naturally solves the coincidence problem. Also, it was found 
that the causal connection scale Rqc consistent with cosmological observations is given 
by Rqq = H + 2H 2 which is proportional to 1Z in a flat universe [21]. This may provide us 
with a physical motivation for the RDE model. Cosmological constraints on this model 
were studied in Ref. [UJ. Similar models have also been studied in Ref. [T2l . 

In this paper, we consider the RDE model with an interaction between dark energy 
and matter, and call it the interacting Ricci dark energy (IRDE) model. We investigate 
the observational constraints on this model obtained from SNIa, CMB and BAO data. 
We organize this paper as follows. In section 2, we describe the IRDE model, and obtain 
analytic expressions for cosmic time evolution. In section 3, we discuss the observational 
constraints on this model. We summarize our results in section 4. 



2 The interacting Ricci Dark Energy model 

We consider the spatially homogeneous and isotropic Universe described by the Friedmann- 
Robert son- Walker metric 

ds * = df 2 _ a 2, t) ( + r 2 M 2 + r 2 gin 2 (2) 

V 1 — kr 2 J 
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where k — 1,0, —1 for closed, flat and open geometries. The time evolution of the scale 
factor a(t) is described by the Friedmann equation 



H = 7^(PA + Pm + P 7 + Pk), (3) 

where pa, Pm, P-y and pu represent energy density of dark energy, matter, radiation and 
curvature, respectively, and H = a/a is the Hubble parameter. 

The energy density of dark energy in the IRDE model is defined as [TD] 



Pa = 3aM 2 (H + 2H 2 + ^), (4) 



where a is a dimensionless parameter. Note that pa is proportional to the Ricci scalar 
curvature 

K = -6 (h + 2H 2 + -^J . (5) 

Moreover, it is assumed that there is an interaction between dark energy and matter. The 
energy densities pa and p m obey the following equations [T%] 

Pa + 3H(1 + w A )p A = -Q, (6) 
p m + 3Hp m = Q. (7) 

The interaction rate is given by [19] 

Q = iHpA, (8) 

where 7 is a dimensionless parameter. The energy density of radiation is given by p 7 = 
p 70 a~ 4 , where p 70 is the present value of radiation density. We adopt a convention that 
a(t ) = 1 for the present age of the Universe t w 14 Gyr. According to eq. (J7|) with Q in 
eq. (jSJ), the interaction can be relevant if 7Pa and p m are comparable, whether or not the 
Universe is in the radiation-dominated epoch. In literatures, various types of interactions 
between dark energy and matter have been considered, including Q = 'yMpH 3 , 'yHpA, 
jHpm, jH(pa + p m ), and so forth. At present, no definite mechanisms to determine the 
interaction are established. Thus, we simply assume eq. ([S]) as a phenomenological model. 
Combining with eqs. (T4]) and ([7]), the Friedmann equation (E]) is written as 



2 dx 2 V 2 2 J dx 



Pl ° -{l- a (l + 1 )}ke- 2x = 0, (9) 



3M2 



where x = In a. The solution to eq. OH]) is obtained as 

H 

r 



= A + e° +x + A.e°- X + A^ x + A k e~ 2x , (10) 



1 If the interaction rate is given by Q — 9b 2 MpH 3 (b 2 > 0) [18 , it turns out that the matter density 
p m becomes negative for a <C 1. This problem does not occur for b 2 < 0. However, the case b 2 < is 
disfavored by observational constraints. 
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where 



2 - la - cry ± J (2 - a) 2 - 2a(a + 2)7 + a 2 7 2 , , 

2a; 

fi 7 o = p-yo/Pco, = —k/Hl and p c o = 3M 2 Hq. Note that cr± can be imaginary for 
sufficiently large a and 7. This implies that there is a parameter region where H 2 has 
oscillatory behavior. However, this region is not phenomenologically viable. The constants 
f2 7 o and Qko are the present value of f2 7 and Qk, respectively. The constants A 1 , Ak and 
A± are given by 

A 1 = fi 70 , (12) 

At = (13) 

«(fT T + 3)fi fc o + 2fi A o-a(l-fi 7 o)(^ T + 4) 

A ± - ± . [14) 

a{<7 + — er_J 

Without the interaction (7 = 0), eq. (fTOjl reduces to the result obtained in the litera- 
ture [in]. In this case, the constants in eq. ( fill are cr + = —4 + 2/a and cr_ = —3. 
Substituting eq. ( TTU1) to eq. P|, the Ricci dark energy density is given by 



Pa 

Likewise, the matter density is 



p c0 J2 a (f + 2 ) A ^ iX - ( 15 ) 



2 

*=+, 



P 

i=+ 



Pc0 {l-«(y + 2)}Ae-. (16) 



The equation of state of dark energy can be found by substituting eq. ( 1151) into the 
following expression: 

„ A = _i_U 7+ J.*iy (17) 



3 \ p\ dx 

In eqs. ( [151) and (TIEi) . the term proportional to e' 7 "^ is dominant in the past a < 1, 
while the term proportional to e u+x is dominant in the future a ^> 1. As an illustration, 
let us consider the case a = 0.45 and 7 = 0.15 which corresponds to a + ~ 0.25 and cr_ ~ 
—3.0. In the past a < 1, the ratio of eq. (|T6|) to eq. ( TT5|) is p m / Pa ~ « _1 (2 + cr_/2) _1 — 1 
w 3.4, while p m /p A w a~ 1 (2 + (J+/2)" 1 - 1 w 0.045 in the future a > 1. 

Notice that both pa and p m include contributions with non-standard time evolution. 
This typically leads to a constant ratio of pa to p m , and it may help to solve the coincidence 
problem. As pointed out in Ref. [10] , the coincidence problem is ameliorated in the original 
RDE model without interaction where pa and p m were comparable with each other in the 
past universe. In this case, pa starts to increase at low redshift, and the ratio p^j p m 
rapidly grows in the future, since p m ~ e~ 3x in the absence of interaction. On the other 
hand, in the IRDE model, the behavior in the past is similar to that in the original RDE 
model, but the ratio p^j ' p m is constant even in the future. 
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3 Observational constraints 



In this section, we study the cosmological constraints on the IRDE model obtained from 
SNIa, CMB and BAO observations. In what follows, we focus only on the spatially flat 
Universe (k = 0). 

The SNIa observations measure the distance modulus \l of a supernova and its redshift 
z. The distance modulus is defined by 

fi = 5l ° gw Wc (18) 
where di{z) is the luminosity distance given by 

d ^ {1 + z) [mry (19) 

We apply the Union data set of 307 SNIa [22], and place limits on the parameters a, 7 
and by minimizing [23] 

X 2 S N = A-^r, (20) 



where 



307 



B = ^ ~ ^ obs (^) (22) 
i=i a tn 

307 1 



C 



<7 

i=l W 



Here, ^ and //obsfe) are the redshift and the distance modulus of the i-th supernova, 
respectively. The corresponding la error is denoted by a^.. 

The CMB shift parameter is one of the least model dependent parameters extracted 
from the CMB data. Since this parameter involves the large redshift behavior (z ~ 1000), 
it gives a complementary bound to the SNIa data (z < 2). The shift parameter R is 
defined as 

R = V^o / -rrhdz, (24) 
Jo H ( z ) 

where z C mb = 1090 is the redshift at recombination, and fl mQ = p m o/p c o is the matter 
fraction at present. We use the value R = 1.710 ± 0.019 obtained from the WMAP5 
data [21]. The CMB constraints are given by minimizing 

ttMB=(^^V, (25) 



0R 



where i? Q bs = 1-710 and or = 0.019. 
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Observations of large-scale galaxy clustering provide the signatures of the baryon acous- 
tic oscillation (BAO). We use the measurement of the BAO peak in the distribution of 
luminous red galaxies (LRGs) observed in SDSS [3]. It gives 



A = 0.469 
where the parameter A is given by 

Ho \ 

H(z BAO ) J 



0.95 



-0.35 



^4 — -\J fi m o 



0.98 



1/3 r 



±0.017, 



Ho_ 

zbao Jo H[z i 



1 2/3 



dz 



and zbao = 0.35. The \ 2 for BAO is 



Xbao 



A- A, 



obs 



(26) 



(27) 



(2f 



where A ohs = 0.469(0.95/0.98)-°- 35 and a A = 0.017 . 



In Fig.l, we plot the probability contours for SNIa (red), CMB (blue), BAO (green) 
observations in the (a,Q\o)-pla,ne in the case without interaction (7 = 0). The la, 2a 
and 3a contours are drawn with solid, dashed and dotted lines, respectively. The joined 
constraints using x 2 — Xsn + Xcmb + Xbao are shown as shaded contours. The best fit 
values with la error are a = 0.36 ± 0.03 and fl\o = 0.68 ± 0.03 with xLm = 326. The 
joint la region is outside the la regions of CMB and BAO. 

The same contour plots as FigJTJin the presence of interaction with 7 = 0.15 are shown 
in FigfSJ Compared to Fig.l, the CMB contours are shifted to the right. The best fit 
values are shifted as a = 0.45 and f^Ao = 0.73. It is remarkable to note that, unlike the 
case 7 = 0, the joined la region is included in all the la regions of SNIa, CMB and BAO. 

Dependence on the interaction parameter 7 is presented in the (7, a)-plane for I^ao = 
0.73 in Fig.3. The best fit values are 7 = 0.15 ±0.03 and a = 0.45 ±0.03 with a + = 0.25, 
a_ = —3.0 and Xmm = 313. It is found that the CMB constraint strongly depends on 
7. As a result, the joined la region appears in the range with sizable 7. This implies 
that existence of interaction between dark energy and matter is favored by cosmological 
constraints in the IRDE model. 



4 Conclusions 

We have considered the IRDE model where the interaction rate is given by eq. OH]). 
We have derived the analytic expressions for the Hubble parameter (fTOjl and the energy 
density of dark energy ( Fl5|) and matter ( fl6l) . Both p\ and p m include contributions with 
non-standard time evolution. We have also investigated current observational constraints 
on this model from SNIa, CMB and BAO observations. In particular, the CMB constraint 
is strongly affected by the interaction. We have shown that a nonvanishing interaction 
rate is favored by the observations, giving a reduction of Xmin- The best fit values are 
tt A0 = 0.73 ± 0.03, a = 0.45 ± 0.03 and 7 = 0.15 ± 0.03. 
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Contour of ^ 2 (RDE,y=0) 
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Figure 1: The probability contours for SNIa (red), CMB (blue), BAO (green) observations 
in the (a, f2Ao)-plane in the case without interactions (7 = 0). The la, 2a and 3a contours 
are drawn with solid, dashed and dotted lines, respectively. The joined constraints using 
X 2 — Xsn + Xcmb + Xbao are shown as shaded contours. The flat Universe {k = 0) is 
assumed in all the figures. 
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Contour of * 2 (RDE,y=0.15) 




Figure 2: The same contour plots as Fig{T]in the presence of interaction with 7 = 0.15. 
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Figure 3: The contours of x m the (j, o)-plane 
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